Robust D-wave Pairing Correlations 
in a Hole-Doped Spin-Fermion Model for Cuprates 
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Pairing correlations are studied numerically in the hole-doped spin-fermion model for cuprates. Sim- 
ulations performed on up to 12 x 12 clusters provide robust indications of D-wave superconductivity 
away from half-filling. The pairing correlations are the strongest in the direction perpendicular to 
the dynamic stripe-like inhomogeneities that appear in the ground state at some densities. An opti- 
mal doping, where the correlations reach a maximum value, was observed at about 25% doping, in 
qualitative agreement with high Tc cuprates' experiments. On the other hand, pairing correlations 
are suppressed by static stripe inhomogeneities. 

PACS numbers: 74.20.De, 74.20.Rp, 74.72.-h 



The nature of high temperature superconductors is an 
important open problem in the area of strongly corre- 
lated electrons jl]. In this context, models for cuprates 
have been extensively used to search for superconduc- 
tivity (SC) arising from a purely electronic mechanism. 
In spite of this effort, the current situation is still con- 
fusing, with few positive reports of ground-state SC in 
electronic models using truly unbiased many-body tech- 
niques. Recently, a simple spin-fermion (SF) model has 
been proposed for cuprates The main advantage of 
the SF-model is its simplicity for numerical studies, while 
still keeping a realistic, unbiased, and non-trivial charac- 
ter. Previous studies have already shown that several of 
its properties, such as magnetic incommensurability, ex- 
istence of a density-of-states pseudogap at the chemical 
potential, and the shape of the Fermi surface, resemble 
experimental data for high Tc cuprates In addition, 
upon hole doping, the ground state exhibits charge in- 
commensurability due to the formation of hole-rich ver- 
tical and horizontal stripes |^ . The goal of this paper is 
to study for the first time the pairing correlations of this 
model, and determine the role that charge stripes play 
in the pairing process. To our surprise, relatively robust 
pairing correlations in the D-wave channel were detected, 
correlated with the presence of dynamical stripes. 

The SF-model is constructed as an interacting system 
of electrons and spins, mimicking phenomenologically the 
coexistence of charge and spin degrees of freedom in the 
cuprates |p|. Its Hamiltonian is given by 



H 



fc? 



h.c. 



J^Si-Si + j'^Si-Sj, 

i (ij) 

(1) 

where Cj^ creates an electron at site i = (ix, iy) with spin 
projection a, Si=^^^ Cj^fTa/jCi/j is the spin of the mo- 
bile electron, the Pauli matrices are denoted by cr. Si is 
the localized spin at site i, (ij) denotes nearest-neighbor 



(NN) lattice sites, t is the NN-hopping amplitude for the 
electrons, J > is an antiferromagnetic (AF) coupling 
between the spins of the mobile and localized degrees of 
freedom, and J' > is a direct AF coupling between the 
localized spins. The density (n) = l — x of itinerant elec- 
trons is controlled by a chemical potential /i. Hereafter 
t = 1 will be used as the unit of energy. J' and J are fixed 
to 0.05 and 2.0 respectively, values shown to be realistic 
in previous investigations The temperature will be 
fixed to a low value: T=0.01, which was shown before to 
lead to the correct high- Tc phenomenology 

To simplify the numerical calculations, avoiding the 
sign problem, the localized spins are assumed to be clas- 
sical (with I Si 1=1). This approximation is not dras- 
tic, and it was already discussed in detail in Ref. 
The model will be analyzed using a Monte Carlo (MC) 
method, details of which can be found in Ref. [|. To 
study the superconducting properties of the system, the 
pair correlation functions, Cu,(r) = (Aj^j.A!"^), are mea- 
sured. The index w indicates D- or extended S-wave 
pairing, and the pairing operator is given by Aj" = 

J2aJ2e=±x±v'^^'^i,c'^i+e,-a, whcrC 6 = 1 for 6 ±X 

and e = ±j/ for S-wave pairing, while e — 1 for e ~ ztx 
and e = — 1 for e = ±y in the case of c?a,2_y2-wave pairing. 

To study the long distance behavior of the pairing cor- 
relations, results on X (A^ = 8 and 12), clusters are 
here presented. These results show that the D-wave cor- 
relations are stronger than S-wave for all the values of 
the parameters studied. A typical comparison between 
the two correlations is shown in Fig. la for {n)—0.75 on a 
12 X 12 cluster The extended S-wave exhibits strong 
fluctuations, while the D-wave results are more robust 
and smoother. This shows that our SF-model captures 
the essence of hole pairing in AF backgrounds, where it is 
well known that d^2_y2 correlations dominate jl]. Thus, 
from now on, we will concentrate only on the behavior of 
the D-wave pairing correlations (DPC). 
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FIG. 1. (a) D-wave (filled circles) and extended S-wave 
(open circles) pairing correlation versus distance r, on a 
12 X 12 lattice at (n)=0.75. (b) D-wave pairing correla- 
tions along the direction perpendicular to the charge inho- 
mogeneities at the densities indicated, (c) The long distance 
behavior of the correlations shown in (b), displayed in greater 
detail, (d) Correlations for x=6, as a function of density, (e) 
Same as (b) but with correlations measured parallel to the 
stripe-like charge inhomogeneities. 

In Fig. lb we present the DPC versus distance (mea- 
sured in the direction perpendicular to the stripe-like in- 
homogeneities) for several electronic densities. At half- 
filling the correlations are very small , i.e. ^ 10~^, at 
the largest distance, but they develop a fairly robust tail 
of order 10~^ at small hole doping (corresponding to an 
order parameter (A;"')^0.1). At these densities expect- 
ing stronger pair correlations would be unrealistic, since 
the low carrier density as well as the previously exten- 
sively documented small quasiparticle weight Z of holes 
in antifcrromagnets suppresses the signal In addi- 
tion, the pairing operator, being nearest-neighbors, is not 
optimized to fit the actual pair size. Ours is among the 
strongest signals for D-wave SC found in unbiased studies 
of realistic high Tc models, and they are even compara- 
ble to those reported in 2-leg ladders 0. The strongest 
correlations are observed at (n) « 0.75, indicating the 
existence of an optimal doping as in real cuprates 
This behavior can be observed in more detail in Fig.lc. 
In Fig. Id the DPC at distance x=6, is shown as a func- 
tion of electronic density, and the existence of an optimal 
doping is again clear. The dip at (n) ^ 0.85, which cor- 



responds to a state that has nearly static stripes, can be 
qualitatively identified to the x=l/8 anomaly observed 
in the cuprates. ||] The correlations in the direction par- 
allel to the stripes are in Fig.le and, surprisingly, they 
are about one order of magnitude smaller than those in 
Fig. lb. This is a remarkable general trend, and the rea- 
sons for this difference are discussed below. 

Previous studies have shown that the system changes 
from an AF insulator to a metal upon doping Q. The 
large suppression of the pairing correlations at half- filling 
is expected due to the absence of holes. However, as 
observed in Fig. lb, an increase of about four orders of 
magnitude occurs when the system is doped and becomes 
metallic. To understand this phenomenon, we will anal- 
ize the properties of the charge and spin configurations 
(MC "snapshots") that contribute the most to the en- 
hancement of the DPC, starting at densities for which 
fairly static stripes are stabilized. 
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FIG. 2. (a) Representative snapshot of the spin and charge 
degrees of freedom on a 12 x 12 cluster at (n)=0.85, in the 
regime of nearly static stripes. The area of the circles is pro- 
portional to the electronic charge, while the length of the 
arrows is proportional to the projection of the localized spins 
on the X-Y plane. When n(i) < (n) the circles are gray, (b) 
Same as (a) but for (n)=0.68. (c) Charge structure factor 
along selected directions in momentum space for the same 
parameters as in Fig. lb. Stripe-like inhomogeneities for all 
dopings are along the x direction, (d) Magnetic structure fac- 
tor along selected directions in momentum space for the same 
parameters as in part (c). 

The spin and charge distribution for a typical MC 
snapshot at (n)=0.85 is shown in Fig. 2a. When the local 
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density is smaller than the average density the circles pro- 
portional to the local charge density are shown in gray. 
As observed, the gray circles determine two fairly static 
horizontal stripes [[lC| . This inhomogeneous charge dis- 
tribution produces a large peak at momentum (0, 7r/3) in 
the charge structure factor N{q), shown in Fig. 2c. It is 
also clear the nearly perfect AF order in the electron rich 
regions (white circles). The magnetic structure factor 
S{q), shown in Fig. 2d, peaks at momentum (tt, 57r/6) in 
agreement with the AF order observed horizontally and 
the incommensurability induced vertically by the stripes, 
which carry a 7r-shift across them. 

The results in Figs. lb, e are averages over MC time, 
and also over all the sites of the lattice. Measurements 
done directly on the individual snapshots were found to 
be similar to the averages, and indicate that the pair 
correlations are weaken along the AF regions. An in- 
termediate value is obtained along the hole-rich stripes, 
but the largest correlation is observed along the direction 
perpendicular to the stripes (see Fig. 3a). 
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FIG. 3. (a) D-wave correlations for the representative snap- 
shot shown in Fig. 2a. The circles indicate correlations along 
row 6 (counting from the bottom of Fig.2a), which corre- 
sponds to an AF array with local charge ~ 0.97. The filled 
squares are correlations along row 3, which has an average 
local charge 0.67. The triangles are correlations in the di- 
rection perpendicular to the stripes starting from row 2 and 
averaged over all the columns, (b) D-wave correlations for 
the snapshot shown in Fig. 2b. The circles indicate corre- 
lations along row 3 (counting from the bottom of Fig. 2b), 
which corresponds to an AF array with local charge ~ 0.78. 
The filled squares are correlations along row 8, which has a 
local charge ~ 0.63, and the triangles are correlations in the 
direction perpendicular to the stripes starting from row 6. 

An analogous effect, but more enhanced, is observed 
for (n)=0.68. Fig. 2b, where four nearly static horizontal 
stripes are stabilized. N{q) (Fig. 2c) peaks at momentum 
(0,27r/3), while S{q) (Fig.2d) at (7r,27r/3). In this case, 
according to Fig. lb and le, the DPC functions are one 
order of magnitude larger in the direction perpendicular 
to the stripes than in the parallel direction. Note also 
that there is a substantial difference in the DPC per- 
pendicular to the stripes corresponding to the densities 



0.68 and 0.85 here analyzed. To shed some light on these 
issues, in Fig. 3b correlations for the snapshot shown in 
Fig. 2b are presented. The correlations (circles) along the 
AF domains parallel to the stripes, with an average lo- 
cal density of 0.78, are very weak. We believe that 
the AF order is responsable for this depletion. The DPC 
are stronger along the hole-richer stripes, as indicated by 
the squares in Fig. 3b. In this case, the local density is 
^ 0.63. However, the strongest correlations, indicated 
by triangles, occur, once again, along the direction per- 
pendicular to the stripes. Along this direction the local 
charge is very inhomogeneous and magnetic incommen- 
surability occurs. These observations lead us to believe 
that local charge homogeneity and AF order do not fa- 
vor D-wave pairing, while local charge inhomogcneity and 
its associated magnetic incommensurability promote it. 
Thus, the difference in the perpendicular correlations for 
(n)=0.85 and 0.68 mentioned above may be related to 
the AF reduction in electron rich regions, as the system 
is doped away from half-filling. 
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(b) <n>=0.75 
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FIG. 4. (a) Representative snapshot of the spin and charge 
degrees of freedom on a 12 x 12 cluster at (n)=0.78, in the 
regime of dynamic stripes, (b) Same as (a) for (n)=0.75. (c) 
Snapshot for the same parameters as in (b) obtained at a 
later time during the simulation, (d) DPC in the direction 
perpendicular to the stripes using two lattice sizes. 

The largest DPC are obtained for densities where 
charge fluctuations appear to be dynamical , according to 
the snapshots observed, and their time evolutions, dur- 
ing several MC runs. This is in agreement with previ- 
ous expectations based on experimental results for the 
high Tc cuprates [ pl| . To further confirm this picture, in 
Fig. 4a a typical MC snapshot for (rt)=0.78 is shown. Al- 
though the charge distribution is clearly inhomogeneous. 
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the electron-rich AF domains are not separated by static 
stripes. We have observed that the shape of these do- 
mains changes substantially during the MC simulation, 
confirming that the charge structures are dynamical. As 
a result, only a small peak at (0, 7r/2) is observed in N{q) 
(Fig. 2c). The magnetic incommensurability, on the other 
hand, exhibits a very sharp peak at (tt, 27r/3) (Fig. 2d). 
As observed in Fig. lb, the pairing correlations are en- 
hanced in the direction parallel to the spin incommen- 
surability, which is mainly vertical based on S{q). Also 
in Fig. lb it can be seen that the pairing correlations are 
the strongest for (n)=0.75. The corresponding snapshots 
are in Figs. 4b and c. As in the previous case, a dynami- 
cal inhomogeneous charge distribution is observed. Rigid 
stripes appear at times during the simulation (Fig. 4b), 
but they become distorted after a few subsequent MC 
iterations (Fig. 4c). A very weak peak at (tt/S, 7r/3) ex- 
ists in N{q) (Fig. 2c), while a sharp peak at (tt, 27r/3) 
(Fig. 2d) is found in its magnetic counterpart. Similar 
behavior was observed at {n)—0.72 and 0.875. 

Our conclusion is that weak charge inhomogeneity and 
magnetic incommensurability appear to be crucial for 
the development of robust DPC, at least within the SF- 
model. Although N{q) seems nearly featureless in the 
regime with the strongest pairing, the states are not ho- 
mogeneous, as observed from the MC snapshots. In fact 
studies performed by us in the non-interacting system 
show that when the charge is uniformly distributed, the 
magnetic structure factor is featureless, and the DPC 
present strong fluctuations and are suppressed fl^ ] . 

A comparison of results on 8 x 8 and 12 x 12 clusters, 
indicates that charge inhomogeneities become more dy- 
namical as the system size increases. For example, at 
(n}=0.75 the stripes appear static on 8 x 8 clusters but, 
as shown, are more dynamical on 12 x 12 ones. In Fig.4d, 
we present the DPC perpendicular to the charge inhomo- 
geneities for (n)=0.75 on an 8 x 8 and a 12 x 12 cluster. 
The figure shows that, on one hand, finite size effects are 
small but, on the other hand, a long range tail starts de- 
veloping on the 12 X 12 cluster and it is not apparent in 
the 8x8 one; its origin appears to be related to the de- 
velopment of fluctuations in the charge distribution that, 
as said above, are observed only in the larger cluster. 

Summarizing, we have found evidence of robust D- 
wave pairing correlations in a doped SF-model, with a 
charge inhomogeneous ground state. Our results indi- 
cate that static AF order and D-wave pairing compete. 
The latter is enhanced when AF is replaced by magnetic 
incommensurability. Static, stripe-like charge inhomo- 
geneities, only decrease the strength of the pairing corre- 
lations, as compared with the effect of dynamic charge in- 
homogeneities. This is reasonable since a "liquid" charge 
distribution should be more favorable to pairing than a 
"crystal" one. In addition, the hole attraction strength 
caused by antiferromagnetism should be maximized near 
half- filling, where two holes form a bound state in, e.g., 
the t-J model As a consequence, an optimal dop- 
ing emerges where AF pairing is still robust, while static 



stripes do not compete with superconductivity. A sim- 
ilar behavior was observed in Ref. |lj], which provide 
evidence that SC can coexist with static stripe order in 
Lai.6-xNdo.4SrxCu04. However, when this coexistence 
is observed, Tc is at a minimum indicating that static 
stripe order competes with SC. According to our results, 
D-wave SC is expected to be maximized when the charge 
inhomogeneities are the most dynamic, a situation that 
appears to occur in doped LSCO. In this situation, mag- 
netic incommensurability manifest clearly as a peak in 
the magnetic structure factor, while N{q) is almost fea- 
tureless. This is similar to the behavior observed with 
neutron scattering in high Tc cuprates [ pl| . 
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